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ABSTRACT 

We compare sensitive H I data from The H I Nearby Galaxy Survey (THINGS) and deep far UV 
(FUV) data from GALEX in the outer disk of M83. The FUV and H I maps show a stunning spatial 
correlation out to almost 4 optical radii (r25), roughly the extent of our maps. This underscores 
that H I traces the gas reservoir for outer disk star formation and it implies that massive (at least 
low level) star formation proceeds almost everywhere H I is observed. Whereas the average FUV 
intensity decreases steadily with increasing radius before leveling off at ^ 1.7r25, the decline in H I 
surface density is more subtle. Low H I columns (< 2M0 pc~^) contribute most of the mass in the 
outer disk, which is not the case within r25. The time for star formation to consume the available 
H I, inferred from the ratio of H I to FUV intensity, rises with increasing radius before leveling off at 
~ 100 Gyr, i.e., many Hubble times, near ~ 1.7r25. Assuming the relatively short H2 depletion times 
observed in the inner parts of galaxies hold in outer disks, the conversion of H I into bound, molecular 
clouds seems to limit star formation in outer galaxy disks. The long consumption times suggest that 
most of the extended H I observed in M83 will not be consumed by in situ star formation. However, 
even these low star formation rates are enough to expect moderate chemical enrichment in a closed 
outer disk. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: individual (MSB) — radio lines: 
galaxies — stars: formation 



1. INTRODUCTION 

In this letter we compare wide-field images of atomic 
hydrogen (H I) and far ultraviolet (FUV) emission, a 
tracer of recent star formation (SF), in the far outer 
disk of the nearby spiral galaxy M83. Star formation 
in the outer disks of galaxies has been the subject of 
numerous studies, including direct optical o bservations 
of massive SF (e.g.,!F erguson ct al. 1998: Leli eyre fc Rovl 
[20001 : Icimilandre et al.ii2001cide Blok fc Waltedl2003n and 
studies of the star-fo rming interstellar medium (ISM) 
using CO emissi on (iBraine et a l. 2007 1 iGardan et aLl 
[2OOI or dust tZaritskvl ll99lTPopescu fc Tuffs! 1200311 
With its large field-of-view and excellent sensitiv- 
ity, the Galaxy Evolution Explorer (GALEX) revo- 
lutionized this field, revealing widespread, extended 
SF in the far outer reaches (i. e., far beyond the 
optical disks) of many galaxies (iThilker et al.l 120051 
200 ^120091 : IGil de Paz et al.ll2005L [2C)07at IBoissier et al.l 
2007f ). At the same time it has long been known 
that m a ny ga l a,xies host extended H I dist ri butions (e.g. 



Bosnia' '1981'; 'Kamphuis & Briggsl 119921; IB aiaia et atl 
1994; Boomsma et al. 2008; Walter et al. 200g) - incliid^ 
ing M83 (e.g.. lHuchtmeier fc Bohnensteng cl 1981). 

Mainly due to the lack of matched wide-field and sensi- 
tive UV and 21 cm observations, it is still largely unclear 
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how these two extended components - star formation 
traced by FUV emission and atomic hydrogen - relate 
to one another. Although the integrated star formation 
rates (SFRs) in extended UV (XUV) disks represent only 
a small fraction of the total SFR, the H I-SFR connec- 
tion at large radii bears on a number of aspects of galaxy 
structure an d evolution, such as the shapes and edges of 
stellar disks (jPohlen fc T ruiilJo 2006) or c hemical enrich - 
ment gradients across galaxies (Gil de Paz et al.ll2007b( ) . 
The consumption of outer disk H I by SF affects the avail- 
ability of this gas for fueling SF in inner galaxy disks, a 
necessary process f or galaxies to sustain SF over cosmo- 
logical tim es (e.g.. iShlosman . Frank fc Begelman' '1989|; 
[Blitz 1996; iBauermeister et al. 20091. Finally, compar- 
ing UV and H I in the extreme (low-density, often low- 
metallicity) environment of outer galaxy disks can illumi- 
nate the limiting conditions for cloud and star formation 
("star formation threshold"). 

In this letter we combine sensitive, large field-of-view 
H I data fro m THINGS (The H I Nearby Galaxy Survey, 
[Walter eFli,1...2008.') with extremely deep FUV data from 
GALEX to study the relationship between H I (the domi- 
nant component of the ISM at large radii) and FUV emis- 
sion in the nearby spiral galaxy M83. Previous, shallower 
GALEX observations showed that M83 hosts one of the 
most promi nent examples of a XUV disk (Thilker et al.l 
[20051 I2OO7I ). This unique data set allows us to trace 
the distribution of H I and FUV emission across a ~ 50' 
field-of-view, corresponding to almost 4 optical radii r25 
(defined as the 25th B-band mag arcsec~^ isophote). We 
use these data to measure the location (§[3|) and amount 
(§[4|) of FUV emission relative to the H I. In Section §[5] 
we investigate the distribution of H I surface densities in 
different radial regimes across M83. 
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Fig. 1. — GALEX far UV (upper left), near UV (upper right), THINGS Hi (lower left) maps of M83 and a combination of smoothed Hi, 
near, and far UV (lower right). The axes show the offset from the field center at R.A. = 13''37™1.2'^ and Dec. = — 29°52'16". The exposure 
time for the UV data is 13.5 ks, which makes it one of the deepest UV integrations on a nearby galaxy to date. The UV maps reveal a 
rich spiral structure within the star forming disk, a diffuse extended envelope of UV emission and multiple arms, arcs and filaments that 
reach out to almost 4 optical radii (i.e., roughly the edge of the map). All these structures show a stunning degree of correlation with the 
Hi map, indicating (qualitatively) a tight spatial correlation of atomic hydrogen and UV emission far into the outer disk of M83. 
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2. DATA & METHODOLOGY 

We use H I data from the NRACQ VLA obtained as 
part of THINGS to construct a map of H I 21 cm emis- 
sion across a ~ 50' field-of-view (Figure [H lower left 
panel). This map has been corrected for primary beam 
attenuation (the VLA primary beam FWHM is r^ 30'). 
To maximize sensitivity for extended emission, we apply 
the "natural" weighting scheme, which yields a sensitiv- 
ity (Icr RMS) of - 0.8 mJy /beam (or - 0.2MoPC-2) 
and a resolution of ~ 13". M83 is known to contain a 
significant amount of its total H I content in the form 
of diffuse, low column density H I i n its outer disk (e.g., 
iHuchtmeier fc BohnenstengellllQSH ). which our interfer- 

''' The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



ometric observations will not pick up by design. Because 
we are focussing on H I emission associated with (local- 
ized) star formation, the potentially missing flux on very 
large (> 15') scales is not a concern in the current study. 
When quoting H I surface densities we include a factor 
of 1.36 to account for helium and heavier elements. 

The GALEX FUV data are part of our GI program 
GI3-050. Due to two bright stars near M83, the GALEX 
observations were carried out in petal pattern mode al- 
lowing for the fleld-of-view to be centered on the galaxy. 
The GALEX ficld-of-view is slightly larger than the area 
covered by the H I map. We thus restrict our analysis to 
the ^ 50' field covered by both maps. M83 was observed 
for 9 orbits (~ 13.5 ks) yielding a FUV map about three 
times as se nsitive as the GALEX N earby Galaxy Survey 
M83 map ()Gil de Paz et al.l[2007ah . We use a gaussian 
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kernel to degrade the resolution of the FUV map to 13" 
to match that of the H I map. In this map, the (median 
based) la RMS scatter of the noise is about 2 x 10~^ rnJy 
arcsec"^. 

We adopt inclinatio n (24° ) and position angle (225°) 
from iTilanus fc AUenl (Il993l). d istance (D « 4.5 Mpc 



from [Karachentsev et al 

(r25 w 7.7') from LEDA (jPaturel et alllm 



(I2004D and the op tical radius 
For details 



of data processing, the conversion of measured intensities 
to physical units (e.g., FUV intensities into estimated 
SFR surface densities) , assessment of un certainties, etc., 
we refer the reader to iBigiel et al.l (|201Cl[ ). 

3. LOCATION OF FUV RELATIVE TO HI 

Figure [1] shows that H I (lower left panel) and UV emis- 
sion (top panels) are detected far into the outer disk of 
M83 (the edge of the maps corresponds to ^ 3.5 r25). 
It also shows a remarkable spatial coincidence of major 
features found in the maps. Figure [2] shows this even 
more clearly: Plotted are (smoothed) H I contours of 
0.2 Mq pc"2 (~ 2.5 X 10^^ atoms cm~2) on top of the 
GALEX FUV map. Almost every H I feature has a cor- 
responding feature in the FUV map. Even very remote 
structures, such as the tip of the extended, western H I 
arm show FUV emission and thus signs of recent star for- 
mation activity. This close correspondence implies that 
the extended H I disk represents the reservoir from which 
the young stars emitting the FUV emission are forming 
(presumably with an intermediate stage of bound molec- 
ular gas). It also implies, quite surprisingly, that massive 
star formation proceeds - if at a low level - almost ev- 
erywhere H I is observed. 

Figure [3] quantifies the coincidence of H I and FUV|3. 
The y-axis gives the fraction of the total FUV (gray) or 
H I (black) flux found in regions (i.e., along lines-of-sight) 
with H I surface densities lower than the value on the x- 
axis. For example, for an x-axis H I surface density of 
2 M0 pc~^, we identify all area (i.e., lines-of-sight) in the 
maps with Shi < 2 Mq pc~^. We then add up the H I 
(FUV) emission from this area and divide this H I (FUV) 
flux by the H I (FUV) flux of the entire map to get the 
fraction of flux found below the H I surface density of 
2 Mq pc~^. For comparison, we also plot (dotted line) 
the fraction of area as a function of H I surface density 
(e.g., in the above example the fraction of the entire M83 
H I map with Eri < 2 Mq pc"^). 

The key point in Figure [3] is that the FUV curve (gray) 
tracks the H I curve (black) much more closely than it 
follows the area curve (dotted) . The area curve is the ex- 
pected distribution for FUV emission if it were randomly 
distributed across the map. The fact that the H I and 
FUV curves are quite similar thus shows that rather than 
being randomly distributed, the distribution of FUV flux 
seems to follow that of the H I. Also, while most of the 
area in our map is covered by low column (or no) H I, 
the FUV comes mostly from lines-of-sight with relatively 
high Ehi- For example, 90% of the area has H I surface 
densities < 2M0 pc'^^ while only - 25% of the FUV 
flux comes from this area. This behavior is expected if 

** For this plot (and the following) we exclude regions near fore- 
ground stars and background galaxies identified from the GALEX 
pipeline source catalogs and from comparing UV and optical imag- 
ing 



FUV is spatially correlated with H I emission and implies 
that a substantial amount of H I (high surface densities) 
seems to be a necessary prerequisite to star formation. 

4. RADIAL TRENDS 

Having seen that FUV and H I emission are remarkably 
coincident, we now ask how the amount of FUV emission 
compares to the mass of H I along a line-of-sight. The 
H I is concentrated into features (the arms and filaments 
in Figures [1] and [2]) with a relatively low filling fraction, 
and we have just seen that FUV emission is largely coin- 
cident with these features. Our approach here is to look 
at the mean ratio of H I-to-FUV along these features as a 
function of radius, asking "Where there is H I, what is its 
average surface density and what is the associated aver- 
age FUV intensity?" . For this, we consider only emission 
from areas that have H I surface densities (in a smoothed 
30" resolution version of the H I map) > 0.5 Mqpc~^. 

Figure [4] shows the average H I and FUV intensities 
along these features. Differences are immediately appar- 
ent: the FUV profile (gray) decreases over many orders- 
of-magnitude from the inner parts into the outer disk 
before it levels off at about 1.7r25 (note that the FUV 
emission is not corrected for internal extinction, which 
is significant at small radii). The FUV emission at large 
radii is still detected with good sensitivity; the profile is 
always well above the gray dotted line, which shows the 
typical 5(7 sensitivity. This supports the qualitative con- 
clusion from Figure [5] that the extended H I structures 
in M83 seem to coincide with FUV emission, indicating 
recent SF activity out to largest radii. In contrast to the 
FUV profile, the H I surface density (black) shows only a 
moderate change with radius. The typical (azimuthally 
averaged) H I surface density is ~ 2 — 4M0pc~^. 

Assuming FUV emission to arise mostly from young 
massive stars forming at a steady rate, we can convert 
the measured FUV intensities into star formation rate 
surface densities, Esfr (using the same procedure as 
IBigiel et al.l |2010|) . Then the ratio of the two profiles 
Ehi/Ssfr (which is proportional to EhiZ-^fuv) is the 
H I depletion time, the time it takes (present day) SF to 
deplete the current supply of H I. This quantity is shown 
as a function of radius in the right panel of FigureJH (note 
that inside '^ 0.5r25 we do not expect Shi//fuv to trace 
the gas depletion time, as both H2 and extinction become 
important). 

The H I depletion time rises from the optical disk 
r25 out to ~ 1.7r25 and then remains relatively con- 
stant across the far outer disk of M83 at about 100 Gyr, 
i.e., many Hubble times. This very long timescale has 
several implications. First, it is much longer than the 
~ 2 Gyr depletion time of the molecular gas observed in 
the in ner disks of galaxies f Bigiel et al.ll2008l : iLerov et al.l 
l2008f l. Assuming bound, molecular clouds form stars 
with roughly the same efficiency in the inner and outer 
disk, this extremely long H I depletion time suggests that 
the process of building these molecular clouds is the bot- 
tleneck for forming stars at large radii. Combined with 
the rather low H I columns that dominate at large radii, 
the regulating and / or limiting factor for star forma- 
tion in these environments may well be that of driving 
the H I to high enough densities to trigger the atomic- 
to-molecular phase transition and subsequent star for- 
mation. 
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Fig. 2. — Smoothed Hi contours of 0.2 Mq pc ^ overplotted on the GALEX FUV map (compare Figure[T] upper left panel). The figure 
illustrates the close spatial correspondence of Hi and FUV emission along the arms and filaments in the outer disk of M83. 
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Fig. 3. — The distribution of FUV flux (gray), Hi flux (black) and area (dotted line) in our maps as a function of Hi surface density. The 
j/-axis shows the fraction of the total Hi or FUV flux or the fraction of area (relative to the entire map) below a given Hi surface density 
(x-axis). Even though a large fraction of the M83 map is covered by little or no Hi (dotted line), most of the FUV light comes from regions 
with relatively high Hi surface densities (gray line). 



The long H I depletion time also implies that the ex- 
tended H I reservoir in the outer disk is long-lived, or at 
least that in situ star formation is unlikely to consume 
the gas. This means that this gas is in principle available 
as fuel for star formation in the inner disk, if it can be 
efficie ntly transported into the inner part of th e galaxy 
(e.g.. IWong eTaII[200l IVollmer fc Ler"wl[20lol) . Single 
dish observations sensitive to the entire H I distribution 
(compare §[5]) suggest that M83 hosts a very large reser- 
voir of H I beyond t he optical disk - as much as 80% of 
the total H I mass (jHuchtmeier fc Bohnenstengell I1981L 
and references therein) . The total mass of this H I reser- 
voir may well exceed ~ IO^'^Mq, which would in principle 
be enough to sustain present-day SF in the inner disk of 
M83 for 2-3 times longer than is possible from the H2 



alone (jLundgren et al.ll200^ . 

The moderate SF activity at large radii should also 
chemically enrich the H I envelope. If the disk remained 
in its present configuration for ~ 5 Gyr, then the inte- 
grated outer-disk SFR ( ~ O.OSMfDyr"^ beyond rg s) and 
a stellar yield of - 0.06 (jKrumholz fc Dekel|[20l0^ imply 
a mass of heavy elements equal to about 1% of the total 
outer disk H I mass (as determined from our map). This 
is only a coarse estimate, but it shows that even the low 
observed outer disk SFR may enrich the outer disk gas 
to a moderate metallici ty. This fits to o bservations of 
emission from hot dust (iDong et al.ll2008| ) and fairly de- 
cent m etallicities ()Gil de Paz et al.ll2007bl : fBresolin et al.l 
[20091) in the XUV disk of M83. 
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Fig. 4. — Left: Hi (black) and FUV (gray) emission averaged in 6.5" wide azimuthal rings (deprojected radial profiles) in regions of 
significant Hi emission, i.e., along tfie arms in the extended disk of M83 (see text). The error bars show the (Icr) uncertainty in the mean 
in each annulus. The left axis provides the mass surface density scale for the Hi profile, the right axis the intensity scale for the FUV 
profile (we note that in order to convert the FUV intensity scale — in units of mjyarcsec^^ - into SFR surface densities Ssfr — in units 
of MQyr^^kpc"'^ - one needs to multiply /puv by f^ 3.46; compare IBigiel et a l. (2010)). The gray dotted line indicates the typical 5ct 
sensitivity for the (averaged) FUV emission in an (outer disk) annulus. The respective Hi sensitivity is below the lower plot limit. Hi 
and FUV emission show quite distinct radial trends: whereas average Hi surface densities remain relatively constant along the filaments 
in the outer disk, the mean FUV intensity continues to drop before leveling off at 1.7r25. Right: Intensity ratio of Hi and FUV radial 
profiles. Converting FUV intensity into Ssfr, this ratio yields the Hi depletion time (right axis). At large radii, this depletion time remains 
relatively constant at about 100 Gyr, i.e., many Hubble times. 



5. DISTRIBUTION OF HI SURFACE DENSITIES 

We have seen that low level star formation is perva- 
sive in the H I arms and filaments in the outer disk of 
M83. We have also seen that the absolute level of this 
star formation is very low compared to that found in 
galaxy disks. What drives this drop in the fraction of H I 
that eventually ends up in stars? A number of authors 
have suggested that a drop in H I column or volume den- 
sity leads to some critica l prerequisite for cloud forma- 
tion not being met (e.g., lElmegreen fc Parravanol 119941 : 
lSchavd[200l JBush et al.l l2010[). We conclude bv brieflv 
examining how H I surface densities actually change be- 
tween inner and outer disk. 

The left panel in Figure [5] shows the (normalized) 
distribution of Efji for lines-of-sight in different radial 
regimes (see caption for details). The right panel shows 
the same data, but now each line-of-sight is weighted by 
mass, so that the left panel gives the distribution of sur- 
face densities and the right panel shows at what column 
most of the mass is found. 

In the left panel, the histograms from different regimes 
are not particulary distinct though one can see some en- 
hancement in the H I columns found inside r25 (visible 
also as the "bump" in the radial profile in the left panel of 
Figure m between ~ 0.2 and ~ 0.7r25). These differences 
are much easier to see in the mass-weighted histograms 
on the right. While in the optical disk H I surface densi- 
ties between ~ 2 — 8 M0 pc~^ contribute roughly equally 
to the integrated H I mass in this regime, it is clearly the 
lowest H I columns that contribute most to the mass in 
the outer disk (i.e., beyond r25). In fact, there appears to 
be a radial trend in the sense that as one moves inwards 
from the far outer disk, an increasing fraction of the H I 
mass comes from lines-of-sight with relatively high H I 
surface density. 

We saw that in the outer disk of M83, high H I surface 
densities, likely a prerequisite for molecular gas and sub- 
sequent star formation, contribute a decreasing fraction 
to the overall H I mass as one moves to large radii. The 



extended arms seen in Figure [T] thus seem to represent 
merely the "tip of the iceberg" of the extended H I dis- 
tribution in M83. These arms, however, appear to be the 
regions where pervasive star formation can proceed even 
at extreme galactocentrlc radii. 
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Fig. 5. — Normalized histograms of Hi surface density (left panel) - the right panel shows the same distribution but mass-weighted - for 
each line-of-sight in 3 radial regimes; — 1 r25 (gray), 1 — 2r25 (black hashed), 2 — 3r25 (black dotted) and for the entire Hi disk of M83 
(thick black line). Wc bin emission between 0.2 (the sensitivity of the Hi map) and 15 Mq pc~ in equally sized bins. The plots show 
that low Hi columns dominate most of the area at all radii in M83 and that while higher columns contribute most to the mass within the 
optical disk, low column Hi dominates the mass budget in the outer disk. 
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